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Generation of high intensity attopulses is investigated in cylindrical geometry by using 3D particle- 
in-cell plasma simulation code. Due to the rotation symmetric target, a circularly polarized laser 
pulse is considered propagating on the axis of a hollow cone-like target. The large incidence angle and 
constant ponderomotive pressure leads to nano-bunching of relativistic electrons responsible for the 
laser-driven synchrotron emission. A numerical method is developed to find the source and direction 
of the coherent radiation that is responsible for the existence of attopulses. The intensity modulation 
in the harmonic spectrum is well described by the model of coherent synchrotron emission extended 
to the regime of higher order 7 -spikes. The spatial distribution of the higher harmonics resembles a 
spiral shape which gets focused into a small volume behind the target. 


The theory of surface high harmonic generation 
(SHHG) on the boundary of an overdense plasma was 
developed and well established a decade ago [H-Q. Two 
main mechanisms are responsible for the generation of 
higher harmonics in overdense plasma: the coherent wake 
emission and reflection from a relativistically oscillating 
mirror (ROM) Q, which is more efficient at highly rel¬ 
ativistic intensities. A weak transition between them is 
observable at non-relativistic laser intensities by varying 
the plasma scale length Q. In the case of near-critical 
plasma isolated attopulses can be produced by the effect 
of phase-dependent plasma mirror deflection Q. The 
universal spectrum proposed for the ROM spectrum in 
Ref. 0 found to be invalid in some cases of oblique inci¬ 
dence or for plasma boundaries possessing finite density 
scale length [10. Few years later a new mechanism was 
identified providing more intense attopulses via electron 
nanobunching and coherent synchrotron emission (CSE) 
[HI. Generally, the SHHG in the relativistic laser plasma 
interaction can be interpreted as laser-driven synchrotron 
emission [l^ by electrons moving collectively in the laser 
field near the plasma boundary. 

The CSE is analytically described in one dimen¬ 
sion and it is found to be very sensitive on laser and 
plasma parameters, which makes difficult the experi¬ 
mental realization in a controlled way. However, elec¬ 
tron nanobunches are observed in 2D simulations for 
larg e incidence angles (> 70 degrees) and p—polarization 
[la [ 1 ^ . In that case the linearly polarized intense laser 
light can be efficiently converted into attopulses with the 
same polarization. In the present work we study this 
attopulse generation process in 3D in temporal and spa¬ 
tial domains and we show that the measured frequency 
spectrum is well described by numerically evaluated CSE 
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spectra in the regime of higher order gamma-spikes [1^ . 
Instead of a flat foil and linearly polarized pulse we use 
cylindrical targets and circularly polarized (CP) pulse 
in order to ensure conditions for p-polarized incidence 
at each point of the inner target surface. This setup 
automatically results in focusing of the generated radia¬ 
tion, which is a differen t ap proach for achieving higher 
attopulse intensities 0,0| . The generation of ultravio¬ 
let vortices [l8l| has been recently proven using higher or¬ 
der Laguerre-Gaussian modes of the laser pulse. Here we 
show that similar attopulse-vortex (or attospirals) can be 
produced with much higher energy conversion efficiency 
via interaction of CP pulse with overdense plasma. Such 
pulses, carrying orbital angular momentu m [l§| . can be 
used for diagnostics in turbulent plasmas 0|, detection 
of electron vortices 0| or structural exploration of mate¬ 
rials by polarization-dependent absorption spectroscopy 


Modulations in the spectral intensity of emitted radi¬ 
ation has been observed in experiments and simula¬ 
tions 1,0, but in both cases these modulations appear 
as irregular fluctuations, therefore it has never been com¬ 
pared to analytical spectra. Here we present, for the first 
time, a direct comparison between the spectrum of co¬ 
herently emitted radiation and analytical CSE spectrum. 
The most attractive feature of this type of spectrum is the 
constant intensity up to a roll-over frequency, which in¬ 
creases with laser intensity, but it is lower than the value 
predicted by the ROM model 0. Above this frequency 
the spectrum has an exponential decay modulated by pe¬ 
riodic structures. 


The simulations presented in this paper were per¬ 
formed in VORPAL (VSim) 3D. The interaction of a few 
cycle CP pulse with cylinder and cone targets is consid¬ 
ered. The laser pulse propagates along the axis of the tar¬ 
gets, thus the case of cylinder corresponds to tangential, 
while the cone target to grazing incidence. For the cone 
target a =12 degrees half opening angle with respect to 
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the longitudinal axis is considered. The cylinder radius 
is Rt = 0.8 with wall thickness of c? = 0.1 /im, which 
is in the non-transparent regime because the plasma skin 
depth is Ip = aoXLUcr/(T^no) < 50 nm [2^. The total 
length of the target is 2 /rm and its density is no = 28ncr- 
The laser pulse is 20 fs long (10 fs FWHM) with Gaus¬ 
sian temporal profile and with a peak intensity = 
10^°lT/cm^ corresponding to a normalized peak ampli¬ 
tude of oq = (/L[10^®lT/cTO^]A|[^m^]/1.37)^/^ = 6.8. 
The laser wavelength is = 0.8 fim and the focal spot 
size of the Gaussian pulse (FWHM) is 2 /rm. The grid size 
used in the simulation is dx =10 nm in each directions 
(see supplemental material [2^ for detailed simulation 
configuration). 

In order to facilitate a better understanding of the elec¬ 
tron dynamics and to ensure higher resolution, the inter¬ 
action geometry is reduced to 2D and presented in the 
supplemental material [2^. By analyzing the electron 
current and velocity distribution near the plasma sur¬ 
face, the source of coherent radiation could be identified 
in space. The 3D mechanism is the same, but an addi¬ 
tional angle will also define the emission angle (0) of the 
radiation due to the azimuthal motion of electrons. Fig.[T] 
shows three distinct electron populations, each traveling 
on different spiral path: (i) inside of the cylinder pushed 
forward by the laser pressure, near the inner side of the 
plasma wall, (ii) inside the plasma wall and (hi) outside 
of the plasma spreading in radial direction. The emission 
of the observed attopulses is related to the inner electron 
bunches, which have the highest density and energy (see 
pictures with 7 e > 20). In the case of cone target the 
electrons form a more compact continuous spiral, which 
is explained by the better electron bunching observed in 
2D (see supplemental material M)- 
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FIG. 1. Electron density projected on the xy plane for 7 e > 6 
(left) and 7 e > 20 (right) in the case of cylinder target (top) 
and cone target (bottom) at the time instance t = 15 fs. The 
dashed lines show the initial plasma surface. 

The relativistic electrons traveling on a spiral path 



FIG. 2. The isosurface of electromagnetic energy density pro¬ 
duced with cylinder (a) and cone (b) targets is shown at t = 20 
fs for the value u = 3 ■ 10^® J/m® (a) and m = 12 • 10^® J/m® 
(b). The frequencies 5 times larger than the fundamental are 
included. 

inside of the con-like target produce radiation accord¬ 
ing to the basic laws of synchrotron emission. Due to 
the relativistic longitudinal velocity the emitted field ob¬ 
served behind the target at a static position will be 
packed in a very short pulse with about a hundred of 
attoseconds duration. The spiral shape is preserved also 
in the attopulse energy distribution, which is shown in 
Fig. m In this figure the energy density, defined as 
u = eo{Ey + E‘^ + c^By + c^B^)j2, of the filtered electro¬ 
magnetic field extracted from the spatial data at a given 
time is shown. The maximum harmonic number associ¬ 
ated with the plasma oscillations is ~ 5, thus fre¬ 
quencies above this value are used in order to include at¬ 
topulses originated from the electron bunches. For com¬ 
parison, the energy density corresponding to the peak 
amplitude of the laser pulse is ul ~ 3.3 • 10^^ J/m^, 
which is close to the peak value measured in the focus. 

Approximately 0.4 % and 0.9 % of the incident pulse 
energy is converted into attopulses for the cylinder and 
cone target, respectively. This low efficiency arises be¬ 
cause the plasma does not interact with the high intensity 
component of the laser pulse. If only the energy of inter¬ 
acting part of the Gaussian beam (between r = 0.4 /rm 
and 0.8 fim) is considered, then the efficiency is about 
3 %. This suggests this efficiency is achievable with an 
annular laser beam intensity distribution, where the in¬ 
tensity in the center of the beam is zero (see simulation 
setup in HI). The broadening of a divergent spiral can 
already be seen at 25 fs in Fig. [3Ktop). However, the en¬ 
ergy density in the focus point can reach the u > 10^^ 
J/m^ value corresponding to latto ~ 10^° W/cm^ reach¬ 
ing the peak intensity of the incident pulse. 

The focusing effect is clearly visible in Fig. [31 where the 
transversal energy density distribution of high frequency 
fields is shown at different time instances. Averaging 
is made in the x direction, which is the summation of 
energy density in grid cells containing u > 10^^ J/m^ and 
divided by the number of these cells. For the cone target. 
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FIG. 3. Averaged energy density distribution viewed from 
the direction of laser propagation corresponding to cylinder 
(top) and cone (bottom) targets. In averaging only the grid 
cells with u > 10^^ i/w? are inclnded. 


the energy is more concentrated in a small circle near the 
laser axis at t = 20 fs which then laterally expands. This 
focus point is at a; ~ 3.2/rm. These pictures also show 
that the attopulses have a ring structure with isocentric 
arrangement, meaning that the angular distribution of 
the resulting radiation is not uniform. 

The propagation direction of the attospiral defines the 
focusing distance, where the maximum peak intensity is 
expected, thus the details of the emission process has to 
be analyzed. For that we use the temporal evolution of 
the current density near the plasma surface. This method 
is widely used in ID geometry in order to calculate the 
harmonic spectra. In order to transpose this technique 
into 3D, the laser-plasma interaction is assumed to be 
rotational symmetric, which is justified by our alignment 
conditions: the axis of symmetry of the target coincides 
with the propagation axis of the laser pulse. Thus, the 
current density depends on x and t, which is equivalent 
with the phase of the laser pulse. The spatio-temporal 
distribution of the current density components is very 
similar to the one presented in the supplemental material 


The transverse electric field at a given position in the 
laboratory frame is calculated by integrating the current 
density as a function of retarded time, t' = t — (x' — 
x) /c. The distance from the axis r = ^is chosen 
arbitrarily and the mathematical form is 


Ep{r,x = xo,t) = 0.5fioc^ / jpir,x', 

Jo 


( 1 ) 


where xq is a chosen point where the emitted light is ob¬ 
served from and jp = jy + is the transverse current 

density in the yz plane. The measured transverse elec¬ 
tric field from the current density is shown in Fig. |4Ka) 
at different distances from the plasma surface. The am¬ 
plitude of the attopulses emitted near the plasma wall is 
much higher than the peak amplitude of the laser field. 


However, in the far-held the attopulse will have lower 
amplitude because of divergence and incoherence. 

The jx component contains information about the 
emission angle of the EM wave at a given position. Ref. 
Q stated that the surface normal (e„) of the emitting 
electron layer is perpendicular to the emitted electric 
held, e„ x E(x,t) = 0, which, in turn, is parallel to the 

net current carried by the electrons, j = ^jx + jp- From 

this, the emission angle can be expressed as 9{x,t) = 
aicia,Ti{Ex/Ep), where the longitudinal component is cal¬ 
culated as Ex = jx^t. This component propagates in the 
radial direction thus it is different in each point along a 
line parallel with x coordinate. Knowing that the time 
and space data is discrete, the opening angle of the wave 
vector with respect to the x axis can be calculated as: 


e{x, t) = arctan =-(2) 

2^i=[l..n\ 3p\^ Til 

where = iAt, n is the number of time steps and At 
(duration of one time step of simulation) is simplihed. 

The frequency components of the emitted radiation 
can be obtained by applying a Fourier transform on the 
calculated electric held from Eq. [TJ A stable pulse struc¬ 
ture only appears if there is a constant phase difference 
between consecutive frequencies. Ideally, the group de¬ 
lay should be constant over the whole spectrum, but it 
is usually not possible. Extremely short pulse can be co¬ 
herently emitted from an extremely thin electron layer 
(nanobunch) [I^, which consists of relativistic electrons 
travelling nearly at the speed of light in the x direction 
{vx ~ c, 7ex S> 1). Using Eq. [21 it is possible to scan the 
spectrum over a wide range of emission angles. 

Emission angles were calculated for every (x, t) points 
using Eq. then the current densities in Eq. m were 
integrated including the points where the angle is be¬ 
tween 9 and 9 -|- S9, where 69 can be arbitrarily small. 
In this way, the temporal prohle of the emitted pulses 
is obtained for each angle and by applying the Fourier 
transformation, the angular dependent spectrum is ob¬ 
tained which is integrated from w = Swi up to a cut-off 
frequency (filtering) to give the attopulse energy within 
a given angular interval. The resultant angular distri¬ 
bution of the integrated power spectrum (full lines) is 
shown in Fig. HKb). 

The Fourier transformed field data contains also the 
phase of the high harmonics. The degree of coherence, 
the standard deviation of the group delay, is used in order 
to characterize the emitted radiation and is defined as: 

(70 (w) = -yT(<^T^(wyT^T^(w))^, (3) 

where Tg = dtj/dio. It is shown in Fig. Id^b) that the ma¬ 
jority of the coherent radiation is emitted near the plasma 
surface between 0.1 and 0.3 radian, where the measured 
energy density is relatively low. This explains the low 
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energy conversion efficiency earlier observed and means 
that only a small fraction of the radiation is coherent and 
survives far from the surface. Farther from the plasma 
surface {z = 2.45/im), the radiation intensity is much 
lower and it is less coherent (red curves). The peaks in 
the energy distribution (Fig.lHb)) at large angles explain 
the ring structure observed in Fig. |3Ktop). 

The scaling of the spectral intensity with the harmonic 
number is very important in the generation of attopulses. 
Ideally, the amplitude of equally separated frequencies 
should be the same in order to produce high intensity and 
short pulses. The length of the attopulse also depends on 
the spectral interval, which is kept after filtering [2^. In 
Fig.isr a) the Fourier transform of the full radiation (dot¬ 
ted line) is compared to the spectrum emitted between 
0.1 and 0.3 radians (continuous line), where the highest 
degree of coherence is observed. The attopulses shown 
in Fig. [2] and [3] are composed by these higher harmon¬ 
ics, which resemble the spectrum derived in the regime 
ofCSEinRef. [H]. 

The general expression for the spectral intensi^ emit¬ 
ted by a bunch of charged particles is given by [26| : 


Huj) 


dte X {ex J(r, t)) exp[iuj{t — e ■ r/c)] 


(4) 


where e is the direction vector, where the radiation is 
observed from, and r is the position vector of the elec¬ 
tron bunch. Since the radiation emitted in the forward 
direction (x) is of primary interest, the vector product 
in the above expression can be replaced by jp. The ar¬ 
gument of the exponential function can then be written 
as ijj{t — x{t)/c). At the ultra-relativistic limit, the ve¬ 
locity components are defined by the direction of mo¬ 
tion because the absolute velocity is almost constant: 

~ c, where Vx = x{t) and Vp are the lon¬ 
gitudinal and transversal velocity, respectively. In 
we have obtained an expression for Vx, which can be ap¬ 
proximated by a polynomial function: Vx ~ u(l —air^”), 
where v = Ua;( 0 ) is the maximum velocity and t is dimen¬ 
sionless. The polynomial expansion of Vx shows that: 
ai ~ G/oq for oo 1. Now the transversal velocity 
(proportional to jp) is expressed as: Vp = \Jv^ — vj. ~ 
■j2aiVT^. These time dependencies have been used in 
llll| as well for n = 1 {jp crosses zero) and n = 2 {jp 
only touches zero), but higher values have not yet been 
considered. 

In our case higher values of n are also permitted if 
the electron bunch is accelerated fast enough to nearly 
the speed of light and propagates together with the laser 
pulse for a relatively long time. The final expression for 
the harmonic spectrum can be formulated with general¬ 
ized Airy functions [l^, or using the Lienard-Wiechert 
potentials the time dependent electric field can be de¬ 
rived mill, then transformed into the frequency do¬ 
main to obtain the curves presented in Fig. [5^. The 
blue dashed line is evaluated according to the numeri- 




FIG. 4. Electric fields measured at xo = 3 fim using Eq. 
© along different lines parallel with the cylinder wall (a). 
In all cases y = 2 fim. In (b) the angular distribution of 
the radiation for the strongest attopulse is shown including 
only frequencies higher than 5 times the laser frequency. The 
dashed lines show the RMS value of the frequency-dependent 
group delay defined by Eq. ®. The color code is the same 
in both pictures. 



FIG. 5. Frequency spectrum (a) of the strongest pulse shown 
in Fig.HJa) and the corresponding normalized attopulses gen¬ 
erated within the angle of coherent radiation for different laser 
intensities (b). The inset in (b) shows the attopulse produced 
during the previous period of the laser pulse. 


cal integration method described in [T^, where oJrs ~ 
(■\/27)^^"~'’^^'^”(ai)^^^” is a fitting parameter. The red 
dashed line is the Fourier transform of the electric field 
derived in , which shows a better agreement with the 
simulations. 

There is a good agreement for n = 4 and the value of 
OJrs, with oq ~ 4 in the interaction region, suggests that 
the average 7 factor should be around 2 , which is close to 
the average value obtained in the electron bunch. This 
value is confirmed by the second fitted curve (red dashed 
line), where 7 = 2.2 is used. The value of n depends on 
laser parameters and needs to be investigated in future 
studies. In the supplemental material, the spectrum ob¬ 
tained for the cone target is presented. In that case high 
intensity part of the laser also interacts with the plasma 
providing higher 7 of electrons. The numerical fitting 
gives a higher values of n, as at higher intensities, the 
bunch velocity u(l — r^") stays close to c for longer time. 

Fig. [SKb) shows the reconstructed (inverse Fourier 
transform of the extracted coherent spectra) attopulses 
with a duration around 100 attoseconds. Two more sim¬ 
ulations have been performed in order to investigate the 
dependence of attopulse parameters on laser intensity. 
The plasma density is increased to 39ncr when higher 
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intensity is used. With lower intensity we could find co¬ 
herently emitted radiation only in the preceding laser 
cycle, which is shown in the inset and has a higher peak 
value despite of the lower intensity at that time. By look¬ 
ing at the spectrum in Fig. [5^ we can introduce an Wdr 
harmonic number at which the first modulation starts 
and the intensity drops significantly. It is reasonable to 
consider constant spectral intensity up to which 

results in attopulse duration tatto ~ (0.21/a;dr)^L and in¬ 
tensity latto = (wdr/2)^/a;o, where is the intensity of 
the fundamental in the coherent radiation. We have veri¬ 
fied numerically the dependence of ~ SuJrs = (3/2)aQ 
for n » 1, where we assumed 7 ~ ao/2. These scalings 
suggest that the attopulse gets shorter and more intense 
by increasing the laser intensity. The dependence of at¬ 
topulse generation and vacuum propagation on the in¬ 
cidence angle is discussed in the supplemental material 

M- 


In conclusion, a new type of efficient attopulse gen¬ 
eration and an alternative way of their focusing is pre¬ 
sented with energy conversion efficiency as high as few 
percents. The peak intensity of the filtered fields in the 
focus can reach the value of the incident one if a cone 
shape target is used. Practically, a hole in a few microm¬ 
eter thick target is also suitable for this process . The 
spectral intensity measured in these simulations show a 
much weaker decay with the harmonic number than in 
the ROM model, which is very favorable for producing 
intense pulses even on the zeptosecond time scale [^ . 
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